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ABSTRACT 


The  Influence  of  horizontal  wind  shear  layers 
upon  descending*  statically  stable  parachutes  has  been 
studied  analytically  and  experimentally.  The  investigation 
was  concerned  with  horizontal  displacements*  response  velocity* 
response  time*  and  angle  of  attack.  Experiments  were  conducted 
with  two  circular  flat  ribbon*  two  rlngslot*  and  one  ribless 
guide  surface  paraohute  model.  The  results  show  that  the 
analytical  predictions  are  a  satisfactory  first  approximation 
of  the  observed  performance  characteristics. 
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I.  INTRODUCTION 


The  motions  of  a  parachute -load  system  In  a  nearly 
vertical  steady  descent  can  be  determined  relatively  easily 
by  coupling  of  the  known  steady  state  aerodynamic  character¬ 
istics  (Ref  i)  with  the  methods  which  lead  to  the  prediction 
of  the  theory  of  dynamic  stability  behavior  (Ref  2).  When 
descending  from  a  higher  altitude,  the  steady  descent  is 
very  likely  to  be  interrupted  by  horizontal  wind  fields 
known  as  wind  shear  layers .  Upon  entering  such  a  shear 
layer,  the  parachute  system  is  subjected  to  an  initial 
angle  of  attack.  Subsequent  motions  will  be  a  combination 
of  angular  and  lateral  displacements  a3  shown  in  Pig  1. 

The  magnitude  of  these  oscillations  and  displacements  will 
depend,  in  part,  on  canopy  type,  equilibrium  velocity, 
damping,  strength  of  the  wind  field,  its  gradient  with 
respect  to  altitude,  and  apparent  mass  effects. 

In  view  of  these  facts,  a  study  was  conducted 
at  the  University  of  Minnesota  to  determine  analytically 
and  experimentally  the  trajectory  characteristics  of  a 
stable  parachute -load  configuration  under  the  influence 
of  shear  layers  with  constant  and  varying  strengths. 

Experimental  investigations  were  limited  to  five 
parachutes,  two  circular  flat  ribbon  and  two  ringalot  canopy 
types,  each  of  15  and  25  percent  total  porosity,  and  one 
ribless  guide  surface  canopy  with  a  cloth  permeability  of 
120.  Parachutes  were  vertically  injected  into  a  hori¬ 
zontal  flow  wind  tunnel.  In  this  manner,  the  wind  initially 
acted  normal  to  the  longitudinal  axis  of  the  parachute 
canopy.  The  cross  wind  velocities  included  a  steady  wind 
field  of  37  ft/sec  and  a  non-steady  layer  varying  nearly 
linearly  between  5  and  45  ft /sec.  The  parachutes  were 
injected  into  the  constant  flow  field  at  20,  30,  and  40 
ft/sec,  and  they  were  injected  into  the  varying  flow  layer 
at  15  ft/sec. 

Under  the  assumption  of  a  steady  rate  of  descent, 
the  constant  and  the  varying  horizontal  flow  fields  may 
also  be  considered  as  steady  and  non-steady  gusts,  when 
recorded  by  an  observer  positioned  on  the  descending 
parachute -load  system. 

The  analytical  study  of  the  flight  path  of  the 
parachute -load  combination  was  pursued  in  view  of  an  in¬ 
vestigation  concerning  the  wind  response  error  of  rising 
spherical  balloons  (Ref  3).  A  solution  has  been  obtained 
in  the  following  chapters  for  angular  and  horizontal  dis¬ 
placement,  wind  response,  and  response  time  for  the  parachute - 
load  system.  In  a  number  of  cases  a  comparison  between 
calculated  and  observed  behavior  has  been  accomplished. 

The  agreement  is,  in  general,  satisfactory. 
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II.  THBCRBTICAL  ANALYSIS 


A.  Derivation  of  Equations 

When  a  descending  parachute  enters  a  horizontal 
wind  field ,  the  parachute -load  system  will  accelerate  as 
It  passes  through  the  shear  layer.  If  the  parachute  motion 
Is  restricted  to  the  x-s  plane,  as  In  Pig  2,  the  velocity 
may  be  written  In  component  form  as: 

Vp  "  Vxpi  +  Vzpk  .  (1) 

The  horizontal  wind  field  can  be  represented 
by  a  constant  or  a  varying  shear  layer,  as  shown  In  Pig  3. 
m  general,  the  wind  profile  Is  written  with  steady  and 
non-steady  terms  as  follows: 

Vw  -  (V**0  +  Vxw)  1  .  (2) 

Introducing  the  wind  gradient  -  dVXw/dh.  one 
can  write  the  non-steady  wind  term  as  VXW  -  y  (h-ho). 
Whereupon,  the  complete  wind  profile  becomes: 

“  [Vxw0  +Tf(h  -  ho)]  1  .  (3) 

It  Is  also  necessary  to  consider  the  relative 
wind  which  Is  the  difference  between  wind  and  parachute 
velocities,  as  shown  In  Pig  4,  or: 

▼r  -  (Vw  -  V  ,  W 

and  expanding  Eqn  4: 

Vr  -  (V„  -  v,p)  i  -  v„p£ ,  (5) 

from  which  the  magnitude  of  the  relative  wind  follows  as : 

-  i 

|Vr|  -  [  (vw  -  Vxp)2  +  vipj  1  ,  (6) 

and  a  unit  vector  In  the  direction  of  the  relative  wind 
becomes : 


V  -  W|*R|.  (7) 

The  motion  of  a  descending  parachute  Is  governed 
by  Newton's  Second  Law,  written  symbolically  as: 

2>  -  .  (8) 
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CONSTANT  VARYING 


3 


whsrti 


E»  -  sun  of  all  external  foroaa 
»t«  mass  of  paraohute  and  load  plus  included  mass 

5*.  ^  external  f cross  acting  on  the  parachute  *"♦ 
the  aerodynaalo  drag,  gravity,  and  the  apparent  mas  efi«ot*  . 
Hathsmatioally,  these  nay  be  expressed  respectively  set 

Ef  -  B  +  f  +  l>  .  (9) 

To  determine  the  aerodynamic  drag,  it  Is  resumed 
that  the  strength  of  the  wind  shear  layer  Increases  grad¬ 
ually  with  altitude.  Renoe,  dyiamio  effects  can  be  negleottd 
and  the  paraohute  will  align  Itself  with  the  relative  wind 
velocity  veotor,  as  shown  in  rig  4.  The  drag  force  is  then 
written  conventionally  asi 

I-  Oof  8  |VR|a  9,  (10) 

or  introducing  Bqn  6  ami  7  Into  Bqn  10,  it  follows:  , 

*  -  ®D  f  S|VS|[(Tv  -  Vip)  £  -  V„p  k]  .  (11) 

The  weight  force  acting  on  the  system  is: 

F  -  Jgk  .  (12) 

The  apparent  mass  force  arises  from  the  transfer 
of  klnetlo  energy  from  the  parachute  to  the  surrounding 
air  (Refs  4  and  5).  Thin  effect  is  only  noticeable  while 
the  paraohute  is  accelerating  or  decelerating  with  respect 
to  the  air  and  is  written: 

*A  *  m‘  HF*  *  f1-) 

where  m'  is  the  apparent  mass  as  determined  from  potential 
theory  and  experiments  (Ref  6).  Ry  Introducing  Bqn  5, 

Eqn  13  becomes: 

m*  (Vw  -  Vxp)l  -  m'  tq,  k  (14) 

The  equation  of  motion  Is  now  written  using 
Bqns  11,  12,  and  14, 

E*  -  mt(Vxp  1  +  Vzp  k) 

-  <*)£  S|VR|  [(Vw  -  Vxp)  1  -  Vzp  Sc]  (15) 

+  mgk  +  m*  £(VW  -  Vxp)  i  -  Vzp  kj  . 
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op  in  toras  of  tho  horlsontol  and  vortical  component 
directions,  respectively! 


®t*xp  -  00  |  3}Vr|(Vw  -  Vxp)  +  »•  (*w  -  txp)  ;  (16) 

«t**P  -  *  °D  f  S|VR|V,p  +  m«  -  *^ep.  (17) 

It  is  further  assumed  that  the  parachute  descends 
at  a  constant  rate;  hence,  vre  -  0  and  Bqn  17  may  be  solved 
for  the  relative  wind. 


Sine*  Eqn  20  is  *  differential  equation  of  Y»  with  respect 
to  h  only,  a  closed  fora  solution  can  be  found  by  standard 
methods  (Ref  7) .  The  solution  is  simplified  with  the  fol¬ 
lowing  abbreviations. 


Equation  20  can  now  be  written  as: 

dv  +  [  FV  -  Q  -  ZY  ]  dY  -  0  (21) 

Equation  21  is  an  exact  differential  equation  which  is 
solved  as  follows: 

Multiply  Sqn  21  by  eK; 

e™  dV  +  Q**  (PV  -  Q  -  ZY)  dY  -  0  .  (22) 

Letting: 

K  -  e”  (w  -  4  -  ZY)  (23) 

» -  e",  (a*) 

Eqn  22  reduces  to: 

NdV  +  MdY  -  0 .  (25) 


Writing  a  total  differential  equation  as 

dO  -  -  0,  (26) 

where  0  la  any  arbitrary  function,  and  comparing  terms 
in  Bqn  26  with  those  In  Bqn  25,  It  1b  evident  that: 

$-eK,  <«rt 

§y  -  e  W  (FV  -  Q  -  ZY)  .  (28) 

Integrating  Bqn  27  yields: 

a- ve  **  +  *(*),  (29) 

and  differentiating  Bqn  29  with  respect  to  Y: 

§7  -  PVe^  +  f'd).  (30) 

Comparing  Bqns  30  and  28,  one  finds: 

pve^  +  f *(y)  -  e™  fy  -  e**  (q  +  (31) 

Which  reduces  to: 

f'(Y)  -  -  (Q  +  ZY)  ent  (32) 

Integrating  Bqn  32  yields: 

f(Y)  -  (-  f*  +  fs)  e**  “  f  e”  +  Cl  ,  (33) 

where  is  an  arbitrary  constant  of  Integration. 

Now  it  is  possible  to  complete  Bqn  29,  since 
f(Y)  is  known.  Hence, 

o  .  veK  +ere  (-  f  if  +  fs  -  f)  +  oi  .  (34) 
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But  0  and  are  both  arbitrary  and  can  be  combined  Into 
a  single  constant,  C2.  Thus,  Bqn  34  becomes: 

0  -  veW  +  e**  (-  f*  +  ^  -  f)  +  C?  ,  (35) 

and  after  rearranging  terms,  Bqn  35  bebomes: 

v  -  +  fy  -  c2  e-py.  (36) 

The  solution  is  completed  by  Introducing  the  initial  para¬ 
chute  velocity  at  the  boundary  of  the  wind  shear  layer  as: 


■  V  at  h 
*P  xp0 


lo» 


that  Is, 


V  -  Vxpo  at  Y  -  0. 


Therefore, 

Cg  ■  Q  f  Z//P  "  vxpQ* 

and 

v  =  (1  -  e  ~**)  +  f  *  +  vxPo  e -PY.  (37) 

After  substituting  the  respective  values  for  Z,  P,  Q,  Y, 
and  V,  Bqn  37  yields  the  horizontal  parachute  velocity  as 
a  function  01'  altitude. 


'xp 


xwQ 


+  If  (h  -  hQ)  +  (V. 


xp0 


“  VXWo^ 


-ft-  (h 

Vzp 


r  arrirr^-^-ho  n 

U  -  e  vzp 


(38) 


Equation  38  can  be  rearranged  as  follows  to  yield  the  wind 
response,  sometimes  called  the  wind  response  error. 
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Vxw0  +  S'  (h  -  h0)  -  Vxp  =  (VrWo  -  VxPo)  • 


e  (h  '  ho)l  +  t  IT  "F 

c  vzp  J 

[i_  eSt-^nr±(h-ho>l« 


or  by  Eqn  3, 

Vw  -  VXD  -  <V: 


vw  -  Vxp  -  <Vxw0  -  ykPo)  f  „  iffirr  J"  (h 
2  L  C  Vzp 

+  /  —  —  JL _S_  (h  -  hn) 

*  m  8  1  -  e  1  o) 


*•>] 

(39) 


Physically  the  wind  response  is  the  velocity 
difference  between  the  parachute  and  surrounding  air.  As 
the  parachute  enters  a  wind  shear  layer,  the  effect  of  the 
velocity  difference  is  to  accelerate  the  parachute.  For  a 
steady  wind,  the  parachute  accelerates  until  it  moves  with 
the  wind,  as  shown  in  Fig'  5-  When  the  wind  field  is  a  varying 
shear  layer,  as  in  Fig  6,  the  parachute  will  accelerate  to  a 
maximum  response,  which  for  large  values  of  h  -  ho  becomes: 


(yw 


V  ) 
xp'max 


(40) 


It  Is  evident  from  Eqn  39  that  the  wind  response 
depends  strongly  on  the  rate  of  descent,  wind  gradient  and 
velocity  differential  between  parachute  and  wind  as  the 
parachute  enters  the  shear  layer.  TTie  apparent  mass  Influences 
the  response  error.  However,  this  effect  may  be  small  for 
highly  porous  canopies ,  (Ref  6) .  In  the  calculations  of 
Section  III,  the  apparent  and  included  masses  were  neglected 
because  of  Insignificance. 


if  the  analysis  were  applied  to  parachutes  with 
specifically  high  apparent  mass  or  meteorological  balloons, 
the  apparent  mass  should  not  be  neglected,  because  the 
apparent  mass  effect  Is  a  driving  force  in  addition  to 
aerodynamic  drag,  thus  effectively  reducing  the  magnitude 
of  the  wind  response. 


Knowing  the  horizontal  wind,  the  rate  of  descent, 
and  the  wind  response,  one  can  also  compute  the  parachute 
inclination,  angle,  ©  ,  which,  theoretically.  Is  eqqal  to 
the  angle  of  relative  wind  as  defined  in  Fig  4.  Mathema¬ 
tically,  the  Inclination  angle  is  expressed  by: 
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ALTITUDE  (h) 


BOUNDARY  OF  SHEAR  LAYER 


This  relation  follows  if  the  parachute  aligns  itself  in¬ 
stantaneously  with  the  direction  of  the  relative  wind. 

The  horizontal  parachute  velocity,  as  expressed 
in  Eqn  38,  can  be  integrated  again  yielding  the  horizontal 
displacement  or  finally,  the  trajectory.  This  is  accomplished 
after  introducing  the  following  terms: 

dX  cLX  ..  dX 
vxp  "  S&  "  S#  ft  "  ~3h  vzp 

Y  -  h  -  hQ 

dY  »  dh 


JSL 


“t  ~  '  zp 
Equation  38  is  now  written  sb: 


«p 


L  zp 


*Po  0-PY  +  /  Jew o 


zp  \  zp 


»t  V: 


g 


dY(42) 


where  Xp  is  the  parachute  horizontal  position.  Integrating 
Eqn  42  leaves: 


if  v2  Vxpo  -PY 

¥T~  Y  ■  w~© 

zp  zp 


The  constant  of  Integration,  C ,,  is  determined 
from  the  Initial  conditions:  3 

Xp  -  0  at  h  -  hQ  .  Y  -  0. 


(43) 


1 

¥ 


(V; 


XpQ 


-  v. 


mt  V! 


XWq 


Therefore: 


Ti*®- trajectory  equation  la  now  written  completely  as: 


V  Q  Vxw.  m.  V 

*n  "  W~  <h  *  h0)  +  V”2  “  ^  "g^  <h  "  ho* 

*  IP  SB  * 


[nL.  V2  m,.  +  m'  V  "1 

<v*p0  -  +  f  ^  -f)  (-S— )  (-§* )J 

r  5^TT  dh  <h  -  ho>l 

U-e^  Tsp 


r  srfv  J-  (.h  -  ho>l 

U-e^  Tsp 

A  superficial  review  of  the  wind  response  and  trajectory 
equations  might  Indicate  that  they  are  Independent  of  canopy 
size,  drag  coefficient,  weight,  and  air  density.  However, 
these  terms  are  all  incorporated  within  the  terminal  descent 
velocity  Vzp»  expressed  conventionally  as  equilibrium  speed 


The  wind  response  can  also  be  expressed  as  a 
function  of  time  by  assuming  that  the  rate  of  descent  remains 
constant  throughout  the  wind  shear  layer.  Thus,  the  time 
required  for  the  parachute  to  descend  from  ho  to  h  becomes: 


Introducing  Bqn  45  into  Bqn  39  and  after  rearranging 
terms,  the  wind  response  takes  the  form: 


It  Is  convenient  to  Introduce  a  nondimena lonal 
time  T  »  t/tr,  where  the  response  time,  tp,  is  the  time 
required  for  the  parachute -load  system  to  approach,  to  a 
certain  degree,  the  steady  state  conditions  after  entering 
the  shear  layer.  In  this  view,  the  response  time  for  & 
constant  shear  layer.  Fig  7,  has  been  determined  for  the 
condition  that  the  wind  response  diminishes  to  5  percent 


16 


\—£.^^y~--TYU0 


VELOCITY.  (V^.Vi,) 


Altitude  or  Time  Versus  Parachute  and  Wind  Velocity 
for  Constant  Shear  Layer 


TIME 


of  the  horizontal  wind  velocity;  or  numerically. 


vxw0  ~  vxp 


0.05  . 


A  constant  wind  shear  layer  Implies  that  the  slope  Jf 
and  assuming  Vxp0  -  0,  Bqn  47  reduces  to: 


mt  + 


m' 


m 


!*£ 

g 


in 


-  v. 


xp 


'XWr 


). 


0, 


(48) 


Subsequently,  the  response  time  Is  determined 
by  Introducing  the  5  percent  wind  response  requirement;  hence. 


w*  +  n*  V  a 

— -f  & 


(49a) 


or 


ra.  +  m1 


m 


g 


(49b) 


Thus,  Eqn  49b  is  the  response  time  for  a  parachute - 
load  system  in  a  constant  wind  shear  layer,  provided  that 
the  rate  of  d»  scent  remains  constant  throughout  the  entire 
layer.  Continuing  the  analysis,  the  nondimens lonal  time 
is  now  found  by  dividing  Eqn  49b  Into  Eqn  48,  yielding: 


T 


t 

T~  “ 


in  E). 

vxwQ 


(50) 


It  is  Interesting  to  note  that  the  nondimens ional 
time  is  a  function  of  the  wind  and  parachute  velocities 
only.  However,  one  mrsl*.  realize  that  the  expression  for 
horizontal  parachute  .ocity,  Eqn  38,  contains  all  of  the 
parameters  of  any  reaj.  problem. 


For  a  gradient  shear  layer,  the  response  time  has 
been  defined  as  the  time  for  the  wind  response  to  reach 
95  percent  of  the  maximum  wind  response  value,  as  shown  In 
Fig  8.  The  analysis  proceeds  by  dividing  the  argument  of 
Eqn  47  by  the  maximum  wind  response,  (V*  -  Vxp)max>  yielding. 


ra  g 


' v.  -  w 


VXWo  ~  VXPo 


vc'-v-J 


xp'max 


(V  rr — j — 

'  w  xp'max, 


(51) 
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Recalling  Eqn  40  and  again  assuming  VXB  ■  0  and 
VXWo  ■  0,  Sqn  51  simplifies  to: 


m.+m '  VBR  . 

t  — s—  -P  A> 

m  g 


^  "  ^xP^pax  ~  1 

-  1 


And  introducing  the  95  percent  wind  response  requirement 
yields  the  response  time  In  a  gradient  chear  layer  to 

fer  - - m - g  (  ^~T“)  - 


(53a) 


or 


3 


ro^-hi’ ' 


(53b) 


The  nondime ns ional  time  for  gradient  shear  layer 
is  obtained  by  dividing  Eqn  53b  into  Eqn  52,  leaving: 


In  view  of  Eqn  49b  and  53b,  it  is  evident  that  the 
response  times  in  steady  and  non-steady  Bhear  layers  are 
identical  to  each  other  and  independent  of  the  strength  of 
the  she&r  layer.  It  is  also  interesting  to  note  that  the 
response  time  increases  when  apparent  masB  terms  are  intro¬ 
duced  as  shown  in  Pig  9.  Disregarding  the  instant  of  the 
first  impact  of  the  parachute  in  the  shear  layer,  the  effect 
of  apparent  mass  is  a  driving  force  for  non-Bteady  winds  and 
a  retarding  force  for  steady  winds;  consequently,  the  response 
time  increases  In  each  case.  For  h  -  ho  sufficiently  large, 
the  relative  accelerations  between  parachute  and  wind  vanish. 
Hence,  the  apparent  mass  force  also  vanishes  and  the  final 
wind  response  becomes  independent  of  apparent  mass. 
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Calculations  have  been  made  for  a  parachute 
entering  a  wind  Bhear  layer  using  the  results  of  Eqns  38, 

4l,  44,  and  49b.  The  computations  arc  based  on  a  wind 
shear  layer  typical  of  those  found  in  the  atmosphere  (Ref  9) 
which  has  a  gradient  Y  ■  0.1  and  a  maximum  velocity  of 
40  ft/sec.  A  shear  layer  800  ft  deep,  shown  In  Pig  10, 
was  used  together  with  parachute  descent  velocities  of 
25  and  50  ft/sec.  The  results  of  these  calculations  are 
cross  plotted  in  Pigs  10,  11,  and  12. 


Calculations  have  also  been  made  for  a  shear 
layer  which  has  a  gradient  Jf  «  0.1  and  a  depth  of  400  ft. 
These  computations  are  shown  for  altitude  vei3us  horizontal 
velocity  in  Pig  13,  with  terminal  velocity  conditions  of  25, 
50,  75,  and  100  ft/Bec.  A  set  of  calculations  has  also 
been  completed  for  Bteady  wind  shear  layers  of  10,  20,  and 
40  ft/sec,  shown  in  PlgB  14,  15,  and  1 6,  respectively. 

In  each  case  the  terminal  velocity  conditions  were  varied 
from  25  to  100  ft/sec. 


Prom  the  development  of  the  preceding  equations, 
it  is  evident  that  the  particular  effects  of  the  free  or 
damped  oscillations  of  the  parachute,  which  obviously  are 
Initiated  when  the  parachute  enters  the  shear  layer,  are 
neglected.  In  view  of  the  fact  that  within  the  scope  of 
this  study,  merely  statically  stable  parachute  systems 
are  considered,  and  remembering  that  in  accordance  with 
Ref  1  the  oscillations  of  such  systems  diminish  very  quickly 
the  presented  method  appears  to  be  acceptable.  It  may  also 
be  stated  that  Ref  1  shows  that  the  horizontal  displacement 
of  a  statically  stable  system  under  damped  oscillations 
Is  small.  Therefore,  the  presented  trajectory  analysis  may 
also  be  considered  a  good  approximation,  when  considered 
from  this  point  of  view. 


III.  EXPERIMENTS 


A.  Introduction 

Since  the  cross  wind  effects  are  of  great  prac¬ 
tical  importance,  laboratory  tests  were* conducted  tm  check 
the  validity  of  the  analytical  results.  For  this  purpose, 
circular  flat  ribbon  and  ringslot  models  each  of  15  and  25 
percent  total  porosity  and,  for  comparison,  a  ribless  guide 
surface  model  with  a  nominal  porosity  of  120,  were  injected 


Pig  10.  Horizontal  Velocity  as  a  Pune t Ion  of  Altitude  of  Point - 
Mass  Parachute  System  at  Various  Terminal  Velocity 
Conditions,  Wind  Gradient  (#  )  =  0.1  ft/sec/ft  (0<h<400), 
Wind  Velocity  (Vw)  =  40  ft/sec  (400<h<800) 
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ALTITUDE  (h)  100  ft. 


Pig  11.  Trajectory  for  Point-Mass  Parachute  System  at  Various 
Terminal  Velocity  Conditions ,  Wind  Gradient  ( )  = 

0.1  ft/sec/ft  (0<h<400),  Wind  Velocity  (Vw)  =  40  ft/sec 
(400 <h<  800) 
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2  4  6  8 

ANGLE  OF  RELATIVE  WIND,  0  (cleg.) 
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Angle  of  Relative  Wind  as  a  Function  of  Altitude  for 
Point -Mass  Parachute  System  at  Various  Terminal  Velocity 
Conditions,  Wind  Gradient  ( #  )  »  0.1  ft/Bec/ft  (0<h<400) 
Wind  Velocity  (Vw)  =  40  ft/sec  (400<h<800) 


ALTITUDE  (h)  100  ft. 


HORIZONTAL  VELOCITY  ,  Vxp  (ft /sec.) 


Pig  13.  Horizontal  Velocity  as  a  Function  of  Altitude  for  a 

Point-Mass  Parachute  System  at  Various  Terminal  Velocity 
Conditions,  Wind  Gradient  ( u' )  “  0.1  ft/sec/ft 
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ALTITUDE  (h)  100  ft 
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ALTITUDE  (h)  100  ft. 


Fig  15.  Horizontal  Velocity  aa  a  Function  of  Altitude  for  a 

Point-Mass  Parachute  System  at  Various  Terminal  Velocity 
Conditions,  Wind  Velocity  (V*)  -  20  ft/sec 
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Into  the  winds t re am  of  a  horizontal  open  test  section  wind 
tunnel.  A  description  of  the  models ,  testing  facilities , 
and  experiments  Is  given  In  Section  IV. 

2he  steady  state  characteristics  of  the  parachutes 
used  are  known  to  be  quite  different  (nef  2).  Hence,  one 
would  expect  corresponding  performance,  differences  in  the 
cross  wind  experiments. 

Comparing  then  the  results  of  the  simplified 
analytical  process  with  those  of  the  experiments,  one  may 
observe  certain  characteristic  differences  which  would  bo 
caused  by  the  damped  oscillations  or  the  dynamic  behavior 
of  the  parachute  systems.  Considering  these  deviations 
in  view  of  the  theory  of  dynamic  stability.  Ref  1,  one 
may  derive  a  refined  analytical  method  based  on  the  analysis 
given  in  this  study,  combined  with  certain  corrections 
extracted  from  the  theory  of  Ref  1  and  from  the  presented 
experimental  results.  The  predictions  of  this  method  would 
be  more  realistic  than  those  available  at  present.  Unfor¬ 
tunately,  certain  limitations  prevented  the  accomplishment 
of  this  phase  of  the  investigation  at  this  tu.me,  but  a 
qualitative  comparison  will  be  presented. 


B.  Experimental  Procedure 

The  fi*e  parachute  types  mentioned  above  have 
been  tested  in  both  constant  and  varying  wind  fields.  A 
steady  wind  tunnel  speed  of  37  ft/sec  was  used  in  combination 
with  parachute  terminal  velocities  of  20,  30,  and  40  ft/sec. 
Experiments  with  terminal  speeds  below  20  ft/sec  were  found 
to  be  inconclusive  because  of  exceptionally  large  angular 
and  horizontal  displacements.  At  speeds  above  40  ft/sec 
the  horizontal  displacement s  were  immeasurably  small  due 
to  the  high  canopy  loadings. 

Varying  or  non -steady  cross  wind  experiments  were 
conducted  with  an  established  velocity  gradient  increasing 
from  5  ft/sec  at  the  top  to  45  ft/sec  at  the  bottom  of  the 
shear  layer.  This  wind  profile  corresponds  to  a  gradient 
of  *  5.  Attempts  were  made  to  inject  models  at  terminal 
velocities  of  15  and  20  ft/sec.  At  20  ft/sec  it  was  found 
that  the  thickness  of  the  varying  shear  layer  was  insufficient 
to  produce  any  measurable  horizontal  displacement,  while 
at  15  ft/sec  the  parachute  oscillations  were  too  large  to 
yield  reliable  results  within  the  observable  range.  Hence, 
experiments  with  the  gradient  shear  layer  were  discontinued. 
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Results 


The  results  of  experiments  In  the  constant  shear 
layer  are  presented  in  Pigs  17  through  28  as  graphs  of  hori¬ 
zontal  displacement,  horizontal  velocity,  angle  of  attack, 
and  angle  of  relative  wind  as  functions  of  the  altitude. 

In  addition  to  the  experimental  curves,  e  curve  showing  the 
analytical  predictions  is  added  to  each  figure.  These 
analytical  curves  have  been  computed  for  the  particular  ex¬ 
perimental  conditions.  The  five  test  models  are  compared  in 
each  figure  except  for  angle  of  attack  versus  time  in  Pig  29, 
which  only  shows  oscillations  of  a  ribless  guide  surface, 
and  15%  porcalty  ringslot  parachutes  which  are  the  most  and 
the  least  stable  parachute,'  respectively. 

In  detail.  Pigs  17,  18,  and  19  represent  the  tra¬ 
jectory  plots  for  all  models  at  terminal  speeds  of  20,  30, 
and  40  ft/sec  in  a  steady  wind  of  37  ft/sec.  It  can  be 
seen  that  all  models  follow  in  principle  the  theoretical 
predictions.  It  should  be  realized,  however,  that  the  com¬ 
parison  is  merely  extended  over  the  initial  part  of  the 
acceleration  process  and  one  may  expect  a  better  agreement 
in  the  later  portion  of  the  trajectory. 

Prom  the  trajectory  curve  and  knowing  the  time 
at  any  parachute  position,  one  can  determine  the  horizontal 
parachute  velocity  as  shown  in  Pigs  20,  21,  and  22.  In 
each  of  the  above  figures  a  respective  theoretical  curve  has 
been  indicated. 

In  the  case  of  altitude  versus  horizontal  dis¬ 
placement  or  velocity,  it  appears  that  the  high  poroBity 
parachutes,  namely,  the  2.5%  porosity  ringslot  and  ribbon, 
agree  more  closely  with  the  theory,  at  least  in  the  upper 
portion  of  the  shear  layer.  This  is  due,  in  part,  to  the 
slow  angular  response  of  these  parachutes  in  a  cross  wind 
and  the  nearly  constant  effective  drag  force  of  these  para¬ 
chutes  under  these  stability  conditions. 

The  angle  of  attack  arises  in  view  of  the  dynamic 
stability  characteristics  of  the  parachutes  in  free  descent 
and  is  simply  the  difference  between  the  parachute  inclination 
single  and  the  angle  of  relative  wind  as  shown  in  Pig  30.  A 
statically  stable  parachute  without  dynamic  effects  or  over¬ 
shoot  would  always  have  a  zero  angle  of  attack.  The  inves¬ 
tigated  parachutes  have  different  aerodynamic  and  dynamic 
characteristics  and  the  Pigs  23  through  28  indicate  different 
angles  of  attack  and  vary  in  their  alignment  characteristics. 
Since  the  rate  of  descent  is  a  consequence  of  the  effective 
drag,  which  in  turn  is  a  function  of  the  angle  of  attack,  the 
rate  of  descent  of  the  various  parachutes  will  vary  during 
the  initial  phase  of  the  parachute  alignment.  Therefore,  the 
experimental  data  presented  in  Pig  23  through  28  should  also 
be  considered  in  view  of  this  aspect. 
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HOR0ONTAL  DISPLACEMENT,  x  (ft.) 


Pig  17.  Trajectory  of  Various  Parachutes, 
(VZp)  -  20  ft/sec.  Wind  Velocity 


Terminal  Velocity 
(Vw)  -  37  ft/sec 


HOR0ONIAL  DSPUACEMENT,  x  (ft.) 


Pig  18.  Trajectory  of  Various  Parachutes ,  Terminal  Velocity  (Vzp) 
«  30  ft/sec.  Wind  Velocity  (V*)  -  37  ft/sec  1 
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Pig  19.  Trajectory  of  Various  Parachutes,  Terminal  Velocity  (VZp) 
-  40  ft/sec.  Wind  Velocity  (Vw)  -  37  ft/sec 
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HORIZONTAL  VELOCITY,  (ftVsec) 


Pig  20.  Horizontal  Velocity  as  a  Function  of  Altitude  for  Various 
Parachute B ,  Terminal  Velocity  (Vzp)  =  20  ft/sec.  Wind 
Velocity  (Vw)  =  37  ft/sec 
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Horizontal  Velocity  as  a  Function  of  Altitude  for  Various 
Iorachutes,  Terminal  Velocity  (Vzp)  =  30  ft/sec.  Wind 
Velocity  (Vw)=  37  ft/sec 
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Pig  22.  Horizontal  Velocity  as  a  Function  of  Altitude  for  Various 
Parachutes,  Terminal  Velocity  (VZn)  =  40  ft/sec.  Wind 
Velocity  (Vw)  «  37  ft/sec 
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!3.  Angle  of  Attack  as  a  Function  of  Altitude  for  Various 
Parachutes,  Terminal  Velocity  (Vzp)  =  20  ft/sec.  Wind 
Velocity  (Vw)  =  37  ft/sec 
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Pig  26.  Angle  of  Relative  Wind  as  a  Function  of  Altitude  for 

Various  Parachutes,  Terminal  Velocity  (VZD)  -  20  ft/sec. 
Wind  Velocity  (V*)  =>  37  ft/sec  P 
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Pig  27.  Angle  of  Relative  Wind  as  a  Function  of  Altitude  for 

Various  Parachutes,  Terminal  Velocity  (VZp)  «  30  ft/sec. 
Wind  Velocity  (Vw)  =  37  ft/sec 
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Pig  28.  Angle  of  Relative  .Wind  as  a  Function  of  Altitude  for 

Various  Parachutes,  Terminal  Velocity  (Vzp)  =  40  ft/sec 
Wind  Velocity  (Vw)  =  37  ft/sec 


Fig  29*  Angle  of  Attack  as  a  Function  of  Tine  for  the  Ribless 
Guide  Surface  and  151*  Kings  lot  Parachute  at  Various 
Terminal  Velocity  Conditions,  Wind  Velocity  (Vw)  -  3?  ft/sec 


In  general ,  the  results  Indicate  that,  as  the 
terminal  veloolty  of  the  paraohuteB  Is  increased,  the 
experimental  observations  agree  more  closely  with  the 
simplified  the  w  as  formulated  In  Section  II.  Figures 
23,  24.  '.'d  Zj  Illustrate  that,  in  general,  the  angle  of 
atta^*'.  cau  *  d  ty  the  parachute  oscillation  is  greater 
at  low  terro.nal  velocities  than  at  high  terminal  velocities. 
The  oscillation  1b  first  exhibited  when  the  canopy  enters 
the  shear  layer.  This  side  gust  causes  an  angular  motion 
which  1b  a  series  of  damped  oscillations.  In  these  exper¬ 
iments  the  observable  process  was  about  one  half  cycle 
because  of  the  Bmall  shear  layer  involved. 

The  restoring  moment  of  each  canopy  will  determine 
the  magnitude  of  these  oscillations  and  the  rate  of  damping. 
Hence,  one  would  expect  a  highly  stable  and  well  damped 
canopy  to  exhibit  less  oscillations  than  other  parachutes. 

As  an  example,  the  ribless  guide  surface  parachute  is 
compared  wth  the  15#  porosity  ringslot  In  Fig  29.  It 
is  apparent  that  the  ribless  guide  surface  dampens  faster 
than  the  ringslot,  as  would  be  expected.  This  effect  Is 
due  to  the  damping  characteristics,  which  are  related  to 
the  stability  derivative  dCm/doc.  Reference  2  Indicates 
that  the  ribless  guide  surface  Is  statically  more  stable 
than  the  ringslot  parachute  and  the  observed  dynamic 
behavior  is  probably  a  consequence  of  the  stronger  sta¬ 
bility  derivative. 

It  is  also  evident  that  in  principle  both  of 
the  highly  stable  parachutes,  namely  the  25#  ribbon  and 
the  guide  surface  parachutes,  approach  In  any  case  the 
theoretical  angle  of  relative  wind  closer  than  the  less 
8 table  parachutes. 

A  further  observation  is  interesting  to  note. 
Namely,  when  the  parachutes  operate  at  a  positive  angle 
of  attack,  they  present  a  larger  drag  area  to  the  hori¬ 
zontal  wind.  Thus,  one  would  expect  that  the  models  would 
be  accelerated  more  rapidly  than  predicted  by  theoretical 
considerations ,  which  operate  with  a  constant  drag  area. 

Ihls  effect  Is  evident  in  the  trajectory  and  velocity  plots. 
Therefore,  the  experiments  have  been  conducted  under  con¬ 
ditions  which  do  not  fully  satisfy  the  assumptions  made  in 
the  analysis,  which  assumed  as  a  basic  condition  that  the 
parachute  aligned  itself  with  the  relative  wind  instan¬ 
taneously  to  a  zero  angle  of  attack. 

The  experiments  show  that,  with  higher  tennlnal 
speeds  of  the  parachutes,  the  initial  angle  of  attack  is 
decreased  and  the  experimental  results  agree  more  closely 
with  theory.  Furthermore,  It  can  be  seen  that  the  most 
stable  parachutes  align  themselves  faster  to  the  resulting 
velocity,  and  the  angle  of  attack  approaches  zero  relatively 
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fast.  Ihls  experience  shows  that  the  aore  stable  parachutes 
follow  in  this  respect  the  theoretical  predictions  best. 

Both  observations  tend  to  support  the  theoretical  findings 
and  ttae  validity  of  the  analytical  approach. 


D.  Conclusions 

The  behavior  of  steadily  descending  parachutes 
under  cross  wind  influence  has  been  analyzed  on  the  basis 
of  a  simplified  equation  of  motion.  A  number  of  predictions 
concerning  horizontal  displacement*  resultant  velocity* 
and  angle  of  deflection  have  been  made. 

Experiments  were  conducted  which  showed*  in  general* 
agreement  with  the  theoretical  predictions  of  horizontal 
displacement  and  velocity.  Particularly  good  agreements 
were  observed  with  the  high  porosity  ribbon  and  rings lot 
parachutes  as  they  initially  enter  the  shear  layer. 

It  was  also  found  that  the  highly  stable  parachutes* 
namely  the  ribless  guide  surface*  high  porosity  ribbon*  and 
ringalot  canopies  align  themselves  fastest  to  the  relative 
wind  and  assume  fastest  the  horizontal  wind  velocity.  In 
practice*  their  oscillations  would  dampen  out  fastest. 

These  statements*  however,  should  be  understood 
to  be  valid  for  the  range  of  the  investigation  performed 
under  this  study.  Speculatively,  it  may  be  Bald  that  in 
the  later  portion  of  the  trajectories,  all  statically 
stable  parachutes  will  align  themselves  closely  to  the 
theoretical  curves. 


IV.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


A.  Introduction 

The  experimental  equipment  needed  for  the  study 
of  cross  wind  effects  posed  some  unusual  requirements*  and 
a  particular  design  with  various  engineering  compromises 
became  necessary. 


B.  Models 

Five  parachute  models  were  used  in  the  cross  wind 
effects  program.  Including  circular  flat  ribbon  and  ringslot 
canopies  each  of  15  and  25  percent  total  geometric  porosity 
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and  a  120  nominal  porosity  ribless  guide  surface  parachute. 
Photographs  and  gore  patterns  for  these  canopies  are  shown 
in  Pigs  31  through  34. 

All  ribbon  and  ringslot  parachutes  were  constructed 
with  a  nominal  diameter  of  16  inches,  while  the  ribless 
guide  surface  had  a  12.6  inch  design  diameter.  The  ribbon 
and  ringslot  parachutes  had  28  goroa  and  suspension  lines 
and  the  ribless  guide  surface  had  only  12  gores  and  suspension 
lines.  Significant  dimensions  of  each  model  are  given  in 
Pigs  32  and  34. 

•Hie  canopy  sizes  were  selected  in  view  of  diffi¬ 
culties  involved  with  the  fabrication  of  small  canopies. 
Generally,  parachute  models  with  nominal  diameters  less 
than  16  Inches  are  difficult  to  construct  and  still  maintain 
reasonable  aerodynamic  characteristics. 


C.  Cross  Wind  Tunnel 


Considering  the  minimum  size  parachute  models 
to  be  employed,  an  open'  Jet  wind  tunnel  with  an  8  ft  high, 

3.5  ft  wide  nozzle  was  designed.  This  Insured  that  the 
parachutes  would  stay  In  the  wind  field  durlfig  their  descent. 
The  8  ft  high  nozzle  would  also  give  the  model  sufficient 
vertical  distance  to  respond  to  the  cross  wind. 

Aerodynamical ly,  the  wind  tunnel  configuration 
was  based  on  experiments  made  with  a  1/12  scale  model 
tunnel  which  particularly  was  used  to  study  means  of  pro¬ 
ducing  gradient  wind  profiles.  Basically,  the  wind  tunnel 
consists  of  three  sections,  namely  the  nozzle,  diffuser, 
and  fan,  as  shown  schematically  in  Pig  35. 

The  nozzle  section  (Pig  30  has  a  contraction 
ratio  of  2.3  to  1.  Horizontal  shelves  at  the  mouth  of  the 
nozzle  serve  to  straighten  the  flow  after  it  passes  ever 
the  turning  vanes.  Using  thea  burning  vanes,  it  is  possible 
to  maintain  a  velocity  gradient  at  the  nozzle  exit.  Figure 
37  shows  how  the  constant  wind  and  gradient  wind  profiles 
were  obtained  using  the  turning  vanes. 

A  diffuser  was  added  ahead  of  the  nozzle  to  steady 
the  flow  after  leaving  the  fan.  This  section  is  essentially 
8  ft  square  and  24  ft  long.  For  simplicity,  the  section 
was  made  with  straight  walls  except  near  the  fan,  where  a 
transition  faring  was  fitted  to  the  fan. 

Hie  wind  tunnel  was  to  provide  a  40  ft/sec  constant 
wind  corresponding  to  a  flow  rate  of  } 0,000  CFM.  In  view 
of  this  requirement,  a  tubular-centrifugal  fan  was  purchased 
from  American  Standard  Corporation  and  coupled  to  a  40  hp 
motor.  Adjustable  inlet  vanes  on  the  fan  allowed  the  flow 
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RIBBON  FAR /CHUTE* 
28  GORE  X=15*/. 


RING5LOT  PARACHUTE 
28  GORE  X-  15°/. 


RIBBON  FARAOHUTE 
28  GORE  A=25V. 


RINGSLOT  PARACHUTE 
28  GORE  *=  25*7. 


Pig  31.  Circular  Plat  Ribbon  and  Ringslot  Parachute  Models 
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Pig  37.  Method  of  Achieving  Various  Velocity  Profiles  in  Cross 
Wind  Tunnel 


rate  to  be  changed  while  the  fan  ran  at  a  constant  speed. 
For  this  purpose,  ser^o  controls  were  arranged  for  conven¬ 
ient  velocity  control. 


D.  Catapult 

The  available  vertical  space  was  insuf x'ieient 
to  attain  the  desired  terminal  velocity  conditions  hv  free 
falling  models.  Therefore,  a  pneumatic  catapult  was  arranged 
capable  of  accelerating  10  lbs  to  100  ft/sec  within  a 
2  ft  stroke.  A  schematic  diagram  of  this  system  1r  shown 
in  Fig  38. 


In  operation,  the  parachute  model  is  attached 
to  the  bottom  of  the  catapult  push  rod,  as  shown  in  Fig  39. 
In  this  manned  the  driving  force  is  transmitted  directly 
to  the  load,  allowing  the  canopy  to  inflate  as  the  system 
accelerates.  Figure  39  shows  the  parachute  ready  for 
launching  with  4  lines  holding  the  canopy  open  for  easy 
air  inflation.  This  method  proved  to  be  very  successful. 

At  launch,  an  electric  device  releases  ths; 
pushrod-plston  assembly  which  then  moves  under  the  force 
from  compressed  air  stored  in  the  reservoir  above  the 
piston.  Bty  varying  the  reservoir  pressure,  one  oan  obtain 
the  desired  acceleration  and  the  respective  tunnel  velocity 
conditions.  Alter  a  stroke  of  approximately  2  ft,  the 
piston  impacts  the  polyurethane  cushion  which  stops  its 
movement  and  absorbs  the  kinetic  energy.  At  this  instant 
the  parachute  and  load  are  released  from  the  pushrod  and 
enter  subsequently  the  wind  shear  layer. 


E.  Experimental  Procedure 

The  primary  objective  of  the  experiments  was  to 
obtain  experimental  data  which  could  be  compared  with  the 
theoretical  predictions.  For  this  purpose  the  parachute 
models  were  injected  into  the  cross  wind  tunnel  at  speeds 
of  20,  30,  tuiti  tO  xo/bbu,  which  were  also  the  equilibrium 
speeds  of  the  system.  These  tests  were  all  conducted  in 
a  constant  wind  layer  of  37  ft/sec. 

As  in  all  experiments,  a  certain  amount  of  error 
is  Introduced  into  each  test.  Hence,  standard  methods  were 
foxlowed  to  minimize  experimental  error  such  as  averaging 
data  and  curve  fitting.  In  view  of  this  procedure,  the 
results  of  Section  III  may  be  considered  to  be  of  acceptable 
accuracy. 


The  testing  arrangement  is  schematically  shown 
in  Fig  40,  while  Fig  41  presents  the  entire  arrangement. 
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PiS  39*  Parachute  Model  In  Launch  Configuration 
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Pig  40.  Experimental  Set-Up  for  Determining  Parachute  Trajectories 
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Fig  4l.  Catapult  Systems  and  Cross  Wind  Tunnel 


®ie  system  is  composed  of  the  catapult,  wind  tunnel,  high 
speed  movie  oamera,  photocells,  and  timing  equipment. 

The  timing  equipment,  shown  in  Pig  42,  consisted  of  a 
light  beam  osolllograph,  frequency  generator,  and  an 
electronic  counter.  The  purpose  of  this  equipment  was 
to  determine  the  lnjeotion  speed  of  the  parachute  as  it 

nonnss  fh*  nhnfnnn  11m 
sr—’-  —  —  — '  *•*'-  « 

Photographic  ooverage  of  the  trajectory  allowed 
one  to  determine  Instantaneous  position,  velocity  and 
orientation  of  the  parachute.  A  sequence  typical  of 
photographs,  shown  in  Pig  43,  indicates  the  rapid  response 
of  the  paraohute -load  system  to  the  cross  wind.  In  general, 
the  canopies  oollapsed  slightly  when  entering  the  shear 
layer,  but  approach  within  a  short  time  the  angle  of  relative 
wind  and  assume  a  more  regular  form.  Since  the  experimental 
wind  shear  layer  was  relatively  thin,  the  parachutes  did 
not  fully  return  to  a  complete  equilibrium,  but  the  shown 
sequence  pictures  are  characteristic  for  the  paraohute 
behavior  in  the  initial  phase  of  the  trajectory. 


Fig  42.  Timing  Equipment  for  DeteraJ 
Velocities 
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